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An evolutionary process of the fast magnetic reconnection in "free space" which is free from any in- 
fluence of outer circumstance has been studied semi-analytically, and a self-similarly expanding solution 
has been obtained. The semi-analytic solution is consistent with the results of our numerical simulations 
performed in our previous paper (see Nitta et al. 2001). This semi-analytic study confirms the existence 
of self-similar growth. On the other hand, the numerical study by time dependent computer simulation 
clarifies the stability of the self-similar growth with respect to any MHD mode. These results confirm 
the stable self-similar evolution of the fast magnetic reconnection system. 

Subject headings: Earth — MHD — Sun: flares — ISM: magnetic fields 



1. INTRODUCTION 

Magnetic reconnection is the most probable 
candidate of the energy conversion process for so- 
lar flares. Such an explosive event in plasma sys- 
tems is very common in astrophysical problems, 
e.g., geomagnetospheric substorms, YSO flares 
and origins of galactic ridge X-ray emission. How- 
ever, entire macroscopic structure of evolutionary 
reconnection is still unclear. 

Actual magnetic reconnection in astrophysical 
systems usually grows in a huge dynamic range 
in its spatial dimension. The initial scale of the 
reconnection system can be defined by the initial 
current sheet thickness. Let us assume it to be 
of the order of the ion Larmor radius (~ 10° [m] 
for solar corona). Eventually, the reconnection 
system will develop into the scale comparable to 
the initial curvature radius of the magnetic field 
lines (~ 10^ [m] for typical magnetic loop of solar 



corona). The dynamic range of the spatial scale 
in evolution is therefore extremely large (~ 10^ 
for solar flares). For the geomagnetospheric sub- 
storms, their dynamic range of growth is also large 
(~ 10^ for substorms). Such a very wide dynamic 
range of growth suggests that the evolution of 
the magnetic reconnection should be treated as 
a development in free space, and individual con- 
ditions of the outer circumstances of evolution do 
not affect the evolution of magnetic reconnection, 
at least in the initial stage. 

The numerical study of magnetic reconnection 
in free space has been performed by our group 
(see Nitta, Tanuma, Shibata and Maezawa 2001), 
and we have discovered a self-similar evolutionary 
process of the fast reconnection. In our work, a 
stable self-similar growth of the reconnection sys- 
tem was shown. We here note the definition of 
the self-similar growth as following: The solution 
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at arbitrary time is the same with the solution at 
any other time if we change the scale of physical 
quantities (see section pl2| ). The spatial scale is 
expanding in proportion to the time. Such self- 
similar evolution is expected to continue unlim- 
itedly in our case. However, simulated dynamic 
range of self-similar growth was restricted (~ 10^) 
owing to the restriction of computer memory and 
CPU time. 

The reconnection model with a pair of slow- 
mode shocks was first presented by Petschek 
(1964). In this model, the magnetic energy 
is mainly converted at the slow-mode shocks. 
Therefore, the time scale of energy conversion is 
determined by that of MHD wave propagation, 
and, as such, is much shorter than that of the 
simple diffusion model (Tajima & Shibata 1997) 
and the Sweet-Parker model (Sweet 1958, Parker 
1963). According to the Petschek model, the 
reconnection rate is almost independent of the 
magnetic Reynolds number (oc log Rm ~ Rrn')- 
Hence, we call this quick magnetic energy conver- 
sion "fast reconnection". 

We should note, however, that the question of 
what controls the reconnection rate is still open. 
Concerning this point, there have been two dif- 
ferent general ways of producing reconnection. 
One assumes that external boundary conditions 
should control the reconnection, so that the resis- 
tivity has no effective influence on the energy con- 
version (see, e.g., Petschek 1964, Priest & Forbes 
1986 for theoretical studies, and Sato & Hayashi 
1979 for numerical simulations). Another way has 
been proposed in a series of numerical works orig- 
inated by Ugai & Tsuda (Ugai & Tsuda 1977, 
1979, Tsuda & Ugai 1977, and recently, Ugai 
1999). In their simulation, they put a localized 
finite resistivity in the current sheet to represent 
an anomalous resistivity that may exist in plas- 
mas. This resistivity acts as a trigger to start 
the magnetic reconnection. Once this resistivity 
is put in, the reconnection starts and evolves self- 
consistently, and forms a fast-reconnection system 
with the Petschek-type structure (a pair of slow- 
mode shocks is formed along the current sheet). 

However, even in the previous numerical studies 
aimed at clarification of the time evolution of fast 
reconnection, the evolution could not be followed 
for a long time. This is mainly owing to the finite 
size of the simulation box, hence the application of 
the results has been limited to spatial scales typ- 
ically, say, hundred times the spatial scale of the 



diffusion region. We are interested in evolution- 
ary process in a free space without any influence 
of the outer circumstance. In such a system, the 
evolution and resultant structure would be quite 
different from these previous models. The key 
process of the self-similar reconnection is summa- 
rized as following. 

We suppose a two-dimensional equilibrium state 
with anti-parallel magnetic field distribution, as 
in the Harris solution. When magnetic diffusion 
takes place in the current sheet by some localized 
resistivity, magnetic reconnection will occur, and 
a pair of reconnection jet is ejected along the cur- 
rent sheet. This causes a decrease of the total 
pressure near the reconnection point. Such in- 
formation propagates outward as the rarefaction 
wave. In a low-/5 plasma (/? -C 1 in the region 
very distant from the current sheet [asymptotic 
region]; as typically encountered in astrophysi- 
cal problems), the propagation speed of the fast- 
magnetosonic wave is isotropic, and is much larger 
than that of other wave modes. Thus, the infor- 
mation about the decreasing total-pressure prop- 
agates almost isotropically fast-mode rar- 
efaction wave (hereafter we call it FRW) with a 
speed almost equal to the Alfven speed Vao in 
the asymptotic region. Hence, the wave front of 
FRW (hereafter we call it FRWF) has a cylin- 
drical shape except near the point where FRWF 
intersect the current sheet. When FRWF suffi- 
ciently expands, the initial thickness of the cur- 
rent sheet becomes negligible comparing with the 
system size V^o^, where t is the time from the on- 
set of reconnection. In such there is only 
one characteristic scale, i.e., the radius of FRWF 
iyAot)-, which increases linearly with proceeding 
of the time. This is just the condition for the 
self-similar growth. 

In this paper, our attention is focused on the 
analytic study of self-similar stage of magnetic re- 
connection in free space. An analytic approach 
has the following importance. 

1) By finding the analytic solution for the self- 
similar growth, we can verify the possibility of the 
self-similar evolution more rigorously. We set the 
initial thickness of the current sheet D. The self- 
similar stage is realized in the limit of Vao^ ^ D. 
However, in the computer simulation, we can only 
perform calculation for a finite duration owing to 
technical reasons. Hence, we will never reach the 
exact self-similar stage. But, if we can verify that 
the result of the numerical simulation is very sim- 
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ilar to the analytic self-similar solution, we can 
accept that the result presented by our numerical 
simulation is truly a self-similar evolution. 

2) Existence of the analytic solution ensures 
that the self-similar growth will continue unlim- 
itedly. This is important because, in the com- 
puter simulation, we cannot continue the calcula- 
tion more than the case we performed in the pre- 
vious paper owing to technical reasons. Of course, 
we should note that when the system sufficiently 
grows to a scale similar to the initially imposed 
system size (e.g., the scale of the flux tube for 
the case of solar flares), this self-similar evolution 
must be modifled by its circumstance. We should 
note that such analytic treatment cannot ensure 
the stability of the self-similar evolution, while it 
is ensured by our numerical simulation. Thus the 
analytic treatment and numerical simulation are 
complementary to each other. 

This paper is organized as follows. In section 2, 
we introduce a special coordinate system called 
"the zoom-out coordinate". In this coordinate, 
the real spatial scale is shrinking linearly as the 
time proceeds. If we choose the origin of time 
appropriately, self-similar expansion can be a sta- 
tionary solution in this coordinate. We should 
modify the MHD equations in a relevant form in 
the zoom-out coordinate. In section 3, we re- 
view basic assumptions which will be made in 
this work. We study magnetic reconnection in the 
limit of very small reconnection rate, and adopt 
the method of perturbative expansion by using 
the reconnection rate small perturbation pa- 
rameter. The initial equilibrium is treated as the 
zeroth order (unperturbed) situation. We solve 
the linearized problem for the flrst order quan- 
tities, which show the variation from the initial 
state. We can rearrange the flrst order equations 
to a single 2nd order partial differential equation 
for magnetic vector potential. This procedure is 
called the Grad-Shafranov approach. We will see 
that this equation is elliptic in our problem. Using 
the SOR routine for solving the Grad-Shafranov 
equation, we obtain a self-similar solution. The 
boundary condition used in this work is discussed 
at the end of section 3. In section 4, properties 
of this self-similar solution are discussed. In sec- 
tion 5, we summarize our results, and discuss the 
properties of this self-similar evolution model. 

2. MHD EQUATIONS IN ZOOM-OUT COORDINATE 
2.1. Basic assumptions 



We focus our attention on the inflow region 
of the reconnection system which is sponta- 
neously evolving without any influence of the 
outer boundary conditions. Macroscopic aspects 
of the evolution are well described by fluid ap- 
proximation, hence, we can use MHD equations 
for this research. In the fast reconnection mod- 
els, the resistivity plays an important role only at 
the diffusion region which is in the current sheet, 
thus, we can neglect the resistive effects in the 
inflow region. In usual astrophysical plasma sys- 
tems, the plasma-/? value is very small (typically 
P ~ 10~^ for solar corona). Hence, we can as- 
sume the system is fllled with non-resistive, and 
pressure-free (cold) plasmas (of course, except in- 
side of the current sheet). 

We study the reconnection starting from an 
equilibrium state of two-dimensional current 
sheet system with anti-parallel magnetic distribu- 
tion, as like in the Harris solution. As discussed 
in the previous section, the self-similar stage is 
realized in a very late period (Vao^ ^ D). This 
means that the initial current sheet thickness D 
becomes negligible in comparison with the sys- 
tem size VAot. We can estimate the only one sys- 
tem size by the fast-mode wave transit scale VAot. 
This is equivalent to assuming an inflnitesimally 
thin current sheet. 

In the case of the Petschek-type reconnection, 
the resultant reconnection rate is the order of 10~^ 
which is consistent with our numerical simula- 
tion for self-similar fast reconnection (Nitta et al. 
2001). We may expect a rather small reconnection 
rate in many cases. This means that variation of 
physical quantities due to reconnection from ini- 
tial equilibrium is very small, and we can treat 
such variation by a perturbation method. 

2.2. Zoom-out coordinates 

"Self-similar solution" may be regarded 
stationary solution in the "zoom-out coordinate" 
which is deflned as 

r' = r/iVAot) , (1) 

where r is the position vector from the recon- 
nection point in conventional coordinate system 
(here after we call this the "flxed coordinate" ) and 
r' is the position vector in the zoom-out coordi- 
nate system. One must note that Vao in equation 
(|I|) represents the Alfven speed at the asymptotic 
region far from the current sheet, and is constant 
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throughout the evolution, r' can be expressed as 
r' — xi + yj + zk, where x-axis is parallel to both 
the current sheet and initial anti-parallel mag- 
netic field, y-axis is perpendicular to the current 
sheet, and 2;-axis is parallel to the current sheet, 
but is perpendicular to the initial magnetic field 
(hence, dz=0 in this two-dimensional problem). 
The origin of this coordinate system is at the re- 
connection point which is chosen to be the center 
of the current sheet (see, section 3.1 of Nitta et 
al. 2001). 

In order to avoid ambiguity of the term "self- 
similar evolution", let us consider the following 
situation. Let the distribution of a physical quan- 
tity, say Q, to be time varying in the fixed coordi- 
nate (i.e., Q = f[r,t]). We measure the distribu- 
tion of Q in the zoom-out coordinate. HQ = g{r') 
is independent of time, we can say that the distri- 
bution of Q is self-similarly expanding at a con- 
stant speed in the fixed coordinate. In our case, Q 
denotes any of the following variables: mass den- 
sity p, gas pressure P, flow velocity v, magnetic 
field B, and any combination of these quantities 
(e.g., magnetic pressure, total pressure, etc.). 

Let us express the MHD variables in a non- 
dimensional form. Note that we give a general 
formalism here, hence gas pressure is remained 
in the following equations (when we solve these 
equations actually, we will take the cold limit and 
neglect the pressure terms). We adopt the follow- 
ing normalization in this paper, 

v^VAo{v' + r') (2) 

P = Po-pV) (3) 
B^VAoVJI^o-B\r') (4) 

P = /3/2.poVAo'-PV) (5) 

where /Iq is the magnetic permeability in vacuum, 
(3 = (Cso/Vao)^ with Cso being the sound speed at 
the asymptotic region. Note that we use SI units 
throughout this paper. 

2.3. Non-dimensional MHD equations in 
zoom-out coordinate 

By using above normalization, the MHD equa- 
tions in the zoom-out coordinate can be put into 
the following non-dimensional form (see appendix 
A), 

V' • {p'v') = -2p' (6) 
V ■ [p'v'v' - B E' + I{/3P' + 5'2)/2] = -3p'^) 



V' ■ {v'B' - B'v') = -B' (8) 
Py-T = 1 (9) 

where / is the unit tensor and 7 is the specific 
heat ratio. 

The left-hand side of each equation is very simi- 
lar to its counter part in ordinary MHD equations, 
but strange source terms appear in the right-hand 
side. These source terms appear as apparent ef- 
fects in the zoom-out coordinate. We should note 
that our main focus is on the structure of the in- 
flow region, which is located outside of the diffu- 
sion region, and upstream of the slow shock as in 
the Petschck model. The resistivity does not play 
any important role in the inflow region. Thus, 
we can use ideal (non-resistive) MHD equations. 
For simplicity, we adopt the polytropic relation 
instead of the full energy equation. The flow does 
not experience a violent entropy production in the 
inflow region (up-stream side of the slow shock). 
Thus, the polytropic variation is a good approxi- 
mation in the inflow region. 

3. EQUATION FOR SELF-SIMILAR SOLUTION 
3.1. Perturbative expansion 

We are treating here a slow energy conversion 
induced by magnetic reconnection. The word 
"slow" used here means very small reconnection 
rate, and should not be confused with that in the 
term "slow reconnection", which usually means 
an explicit dependence of the reconnection rate 
on the magnetic Reynolds number. As we found 
numerically, the self-similar reconnection we wish 
to study is "fast" reconnection in the sense that 
the reconnection rate does not depend upon the 
magnetic Reynolds number (at most, it depends 
logarithmically) . 

In such deviation owing to reconnection 

from the initial state is very small in the inflow 
region. Hence we can adopt the perturbative ex- 
pansion method. We treat the quantities of the 
initial equilibrium state to be the zeroth order, 
and any variation from it should be treated as 
the first order quantity. In this paper, our atten- 
tion is focused to fundamental properties of the 
self-similar reconnection; thus we treat up to the 
first order variation. We can recognize that the 
small expansion parameter in this work is the re- 
connection rate. 

We assume an anti-parallel two dimensional 
magnetic distribution with infinitesimally thin 
current sheet. Note that the self-similar stage 
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is established in the case that there is no fixed 
proper length in the system. However, the sys- 
tem does have a fixed scale: the initial thickness 
D of the current sheet is finite. Thus, such self- 
similar stage is realized only if VAot ^ D (see 
section 2 of Nitta et al. 2001). 

This initial equilibrium state (uniform distribu- 
tion) is expressed as 



B' 



Vo = -r 



P'o 



(10) 

(11) 
(12) 

(13) 



in the upper half plane y > 0. The quantities 
having suffix are the zeroth order quantities. 

The first order quantities represent the devia- 
tion from the initial equilibrium. By substituting 
the expansion form of each quantity 



B' = B'o + B[ 

v' = v'q + v[ 
P' = Po + Pi 
P =Po + Pi 



(14) 
(15) 
(16) 
(17) 



together with equations (|TUp - ( [TB| ) into the MHD 
equations (§l)-(^ in the zoom-out coordinate, we 
obtain the equations for each order of magnitude. 
The zeroth order equations are satisfied automat- 
ically since the zeroth order quantities show the 
uniform equilibrium state (see appendix B). 

In order to solve the reconnection in a low (3 
plasma, let us take here the cold limit /3 ^ 0. 

We obtain the first order equations (see ap- 
pendix B), 



r' ■ vv; = 



(18) 

-r' ■ V'v[ - i ■ VB[ + V{i ■ B\) = (19) 
r' ■ V'B\ + i ■ Wi - r' ■ {V'p[)i = (20) 



Pi = IP'i 



(21) 



Note that even in the cold limit, the dimension- 
less gas pressure P' has a finite value (see equation 
0). When we obtain the value of dimension-less 
density p', we can easily obtain P' by equations 
(|T6|) and (|2l|) . By solving these first order equa- 
tions, we can obtain the self-similar solution. 

3.2. Linearized Grad-Shafranov equation 



We adopt the method of the Grad-Shafranov 
(G-S) equation. Above basic equations (|T8|)-(|2ll) 
show that every physical quantity can be ex- 
pressed function of the z-component A[ of 
magnetic vector potential. Basic physical quan- 
tities are related to A[ in the following way (see 
appendix C). 



Pi 







dx 
dA\ 
dy 

dy 



dy 



■■1- 



(22) 
(23) 

(24) 
(25) 
(26) 

y- 



where the suffix x or y denotes the x- and 
component of vector quantities, respectively. 

By substituting these functional form to the 
above linearized MHD equation (p!8D -(^ID, we can 
obtain the G-S equation for this linearized situa- 
tion (see appendix C), 



:i-x' 



^ - 2xy^^ + (1 



dx'^ 



dxdy 



y 



d^A[ 







dy'^ 

(27) 

In our zoom-out coordinate, the Alfven wave- 
front (or the fast-mode wavefront in the cold 
plasma limit /? — > 0; FRWF) emitted from the re- 
connection point is a unit circle. Hence, we solve 
the G-S equation inside the FRWF. 

We should note that equation (|27|) is a Tri- 
comi type second order partial differential equa- 
tion. We call the region \r'\ < 1 of the upper half 
plane {y > 0) as the region R, and |r'| > 1 as -R. 
The equation is elliptic in R, and hyperbolic in R. 
The region R is affected by the FRW emitted from 
the diffusion region. In R, there is no difference 
from the initial equilibrium state (uniform distri- 
bution), because no signal propagates here yet. 
We do not need to solve the region R. Hence, in 
the following, we solve the elliptic equation (pTl) 
for R under boundary conditions on FRWF and 
on the bottom boundary {y = 0, which is the 
junction surface to the reconnection jet). 

It is convenient to rewrite the G-S equation (pTl) 
and other relations (PB|)-(PB|) in the polar coordi- 
nate (r, 6), because the boundary of FRWF has 
a circular shape. We denote by r = y/x'^ + y'^ the 
distance from the reconnection point x = y = 
and by 6 the angle in the x — y plane from x-axis 
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(i.e., tan 6' = y/x). The G-S equation and other 
relations are rewritten as 



\x 



-— + 1)^0 



2, d'A[ d'A[ dA[ 



dr 



(2^ 



and 





dr ^ 
sin6'^:r^ + — ^ 



dr 



89 



Pi = ip'i ■ 



(29) 
(30) 

(31) 
(32) 



We solve the G-S equation ( Pq ) numerically in the 
region r<l,0<6'<7r with boundary conditions 
on r = 1 and 6' = 0, vr (see next subsection). 

3.3. Boundary conditions 

At FRWF (r = 1), deviation of any quan- 
tity from the initial equilibrium state must van- 



H and [H-il). Thus, 



ish (see equations [[ 

A'-^ = const., and should be a C^-class continuous 
function at r = 1. Without any loss of generality, 
we can set A[ = and drA[ = on r = 1. 

Note that if we impose a Dirichlet type bound- 
ary condition A'^ = on r = 1, another condition 
drA[ = (this means B[ = on r = 1) is au- 
tomatically satisfied. This is easily understood 
from equation (PSI). The first two terms in the 
left-hand side vanish on r = 1, then drA[ must 
vanishes if these quantities are continuous at this 
point. 

The boundary condition at ?/ = is not trivial 
as at r = 1. This boundary corresponds to the 
slow shock, and is the interface between the in- 
flow region and the shocked region (reconnection 
jet). We need precise physical discussion to de- 
termine the junction condition for this boundary. 
Instead of solving this difficult problem, we set the 
boundary condition at this interface by adopting 
the result of our numerical simulation (Nitta et 
al. 2001). 

Our previous simulation gives information for 
every quantity just outside the slow shock. We 
adopt the following simple functional form of A[ 
approximating the boundary values of A'^ from 
our previous simulation. 



(for |x| < Xm) 

— — (|x| - l)Ao 

-L Xm 

(for Xm < \x\ < 1) . 



(33) 



(34) 



A[{\xly = 0) 



where Aq (^ 1) is the value of A[ at the reconnec- 
tion point X = y = (this shows the reconnection 
rate v[y at x = y = 0: see equation [§5)5 is 
the location of the contact discontinuity (the in- 
terface between the original current sheet plasma 
and the reconnected inflow plasma), and Xm is the 
location of the minimum of A[. Thus, we have set 
a complete Dirichlet type boundary problem for 
the linearized G-S equation (P?]). We set Aq, Xc 
and Xm to be consistent with our numerical simu- 
lation. For example, the figures shown in the next 
section are for the case Aq = 0.055, Xc = 0.64, 
Xm = 0.84. 

We compare the model boundary-value defined 
above and the simulation result (Case A of Nitta 
et al. 2001) of A[ along the slow shock in figure |I|. 
The solid and dotted lines show the approximated 
model adopted here and the simulation result, re- 
spectively. The actual slow shock has a small tilt 
angle from the x-axis, so that, when plotting the 
simulation result (dotted fine), we have projected 
the A[ values obtained at the slow shock onto the 
The two straight line-segments constitut- 
ing the model boundary condition have been de- 
termined in such a way that they pass the simu- 
lated values at the X-point (x = 0, A[ = 0.055: 
this value represents the reconnection rate), at the 
contact point {x = 0.64, A[ = 0: the vanishing A[ 
corresponds to the contact point between the ini- 
tial current sheet plasma and the infiow plasma), 
and at the minimum A[ {x = 0.84). We can see 
that the solid line approximates the dotted line 
reasonably well. 



4. SELF-SIMILAR SOLUTIONS 

By using the SOR routine, we numerically solve 
the linearized G-S equation (p8|) under the Dirich- 
let type boundary condition as discussed in the 
previous section. A typical result {Aq = 0.055, 
Xc = 0.64, Xm = 0.84, see the last section) is ob- 
tained in the figures || {A[), || {v[y), and || {p[). 

Let us compare the results of this semi-analytic 
study with the result of our simulation in Nitta et 
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al. 2001. Figures ^, ^ and ^ are the semi- analytic 
results showing the first order quantities solved 
here: magnetic vector potential (2;-component), 
velocity (y-component), and mass density, respec- 
tively. Figures ^ and ^ show the numerical 
result in Nitta et al. 2001: the variation from 
the initial equilibrium of magnetic vector poten- 
tial (z-component), velocity (y-component), and 
mass density, respectively. Each quantity is renor- 
malized in the way explained in section 2.2 (see 
[|l|-[Q). Small numbers affixed to some of the con- 
tours denote the level values at these points. We 
can easily find that the result of the semi-analytic 
study is consistent with the result of our numeri- 
cal simulation not only in the topological sense of 
contours, but also quantitatively. 

We have here succeeded to verify that 1) our 
result of numerical simulation (Nitta et al. 2001) 
truly represents a self-similar growth, 2) this self- 
similar growth is stable over a long duration. The 
spatial dynamic range of the numerical study is 
limited (at most ~ 10^) by technical restriction ( 
e.g., restriction of memory and CPU time), but 
the fact that the result of numerical simulation 
is consistent with the semi-analytic solution sug- 
gests that the self-similar growth will continue 
indefinitely. On the other hand, the stability of 
semi-analytic solution with respect to any MHD 
mode is certified by numerical simulation. We can 
say that our semi-analytic study and numerical 
simulation are complementary to each other. 

5. SUMMARY & DISCUSSION 

5.1. Summary 

Together with our time-dependent numerical 
simulation (see Nitta et al. 2001), the linearized 
perturbation solution discussed in this paper has 
ensured the existence of self-similar growth of fast 
magnetic reconnection. 

Thus we propose here the new model describing 
"self-similar evolution of fast reconnection" . The 
time dependent simulation directly solving the 
MHD equations numerically is effective to check 
the stability of the evolving system. However, the 
duration of simulation is restricted owing to the 
restriction of computer memory and run time, or 
stability of the simulation code itself. Hence, even 
if we find the self-similar (-like) behavior in the 
result of simulation, we cannot be convinced that 
this behavior is a true one that may continue in- 
definitely. 

On the other hand, an analytical study has the 



following properties. If one solves the MHD equa- 
tion in the zoom-out coordinate under the sta- 
tionary assumption and obtain a (semi-) analytic 
solution which is identical to the solution of nu- 
merical simulation, we can be convinced that the 
numerically obtained self-similar-like solution is 
truly the self-similar one. On the other hand, an 
analytical study of the stationary equation does 
not give information on the stability of the solu- 
tion. 

From these arguments, we can reach the fol- 
lowing conclusion. Our time-dependent numeri- 
cal study and analytical study are complimentary 
to each other in establishing the model of "self- 
similar evolution of fast reconnection" . 

5.2. Consistency between semi-analytic and 

numerical studies 

In section 4, we compared our semi-analytic re- 
sult with numerical result. We can notice signifi- 
cant similarity between semi-analytic and numer- 
ical results. 

If we compare them quantitatively in detail, we 
can notice that these two results are consistent in 
the region in which v[y < 0. However, in the re- 
gion in which v'ly > 0, the value of semi-analytic 
result is somewhat different from the numerical 
result. In this region, strong fast-mode compres- 
sion ("piston effect" of the reconnection jet) takes 
place, and, strictly speaking, our linearized treat- 
ment might not be suitable to be applied here. 

Thus, we can conclude that the inflow region is 
well described by our linearized theory. This is 
obviously due to the very small reconnection rate 
(~ 10-2, see Nitta et al. 2001). However, the 
mechanism leading to such a small reconnection 
rate is still unclear. In this self-similar reconnec- 
tion model, the reconnection rate should be self 
consistently determined by the dynamical evolu- 
tion process, and this point will be discussed in 
our forthcoming paper. 

5.3. Condition for self-similar reconnection 

The above type of evolutionary reconnection 
will be realized in the systems in which 1) the 
initial spatial scale of the disturbance ( e.g., scale 
of microscopic instabilities which will lead to the 
anomalous transport phenomena ~ 10° [m] if we 
estimate it by ion Larmor radius in the typical 
case of solar flares) is much smaller than the en- 
tire spatial scale of the system ( e.g., the cur- 
vature radius of the loop of magnetic flux tube 
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~ 10^ — 10^ [m] in the typical case of solar flares), 
and 2) there is no proper spatial scale except the 
one which expands as time proceeds (e.g., the ra- 
dius of FRWF in our case) . In such a system, the 
magnetic reconnection triggered by the initial dis- 
turbance can evolve as the FRW propagates (see 
section 2 of Nitta et al. 2001). The spatial scale of 
FRW propagation is the only proper scale of the 
system if it is much larger than the initial current 
sheet thickness. The spatial dynamic range of the 
evolution is determined by the ratio of the entire 
system scale to the initial disturbance scale. In 
the above mentioned case of typical solar flares, 
it will be of the order of 10^. In such a system, 
the external boundary does not affect the evolu- 
tion for a certain amount of time after the onset 
of reconnection. This means that the system can 
freely evolve independent of the outer boundary 
condition. Such a very wide dynamic range and 
only one evolving spatial scale strongly suggest 
the possibility of self-similar solution. In fact, we 
have obtained the self-similar solution both nu- 
merically and analytically. 

5.4. Near FRWF 

The region near the FRWF (r ~ 1) actually 
has a complicated nature. The region near the 
spearhead of reconnection jet (roughly to say, 
r~l, O<0<7r/8)is influenced not only by 
the fast-mode rarefaction wave, but also by the 
fast-mode compression wave induced by the "pis- 
ton effect" of the jet. Hence one might think that 
the word "FRWF" (fast-mode rarefaction wave 
front) is somewhat misleading. 

However, we should notice that, in our lin- 
earized treatment, any non-linear interaction be- 
tween waves is completely neglected. Thus the 
fast-mode rarefaction wave emitted from the 
vicinity of reconnection point keeps a circular 
shape which is truly located at r = 1 in the zoom- 
out coordinate even if the compressional mode is 
superimposed. From this reason, we may con- 
tinue to call it "FRWF" as we used this term in 
our previous paper. 

5.5. Relation to reconnection jet 

We have treated only the inflow region in this 
paper in spite of the importance of the reconnec- 
tion jet. This is owing to the mathematical diffi- 
culty to treat the reconnection jet where deviation 
from the initial state is very large and the lin- 
ear analysis we have used here breaks down. The 



complete solution which covers the entire region 
has not been obtained yet. 

However, we should note the importance to 
study the property of the inflow region. As Priest 
& Forbes (1986) discussed, the property of the 
inflow region crucially shows the relation to pre- 
vious classical models of reconnection. We can 
see that our new model is an extension of the 
Petschek model to the evolutionary model in free 
space (see section p.7| ). The fast-mode rarefac- 
tion dominated inflow restricts the reconnection 
rate as in the Petschek model (see Priest & Forbes 
1986 or Vasyliunas 1975). In our case, the region 
near the X-point is filled with the fast-mode rar- 
efaction dominated inflow. Hence we may expect 
that our model has a maximum value of the recon- 
nection rate. This point will be clarified by our 
forthcoming study of the reconnection jet which 
induces the fast-mode rarefaction wave. 

In our analytic work, the boundary condition at 
y = 0, which corresponds to the junction condi- 
tion to the reconnection jet, is artificially imposed 
approximating the result of our numerical simula- 
tion. Needless to say, this boundary condition is 
more important than other boundary conditions 
(e.g., conditions on r = 1 or = 7r/2) because it 
represents physical information about the recon- 
nection jet, and it crucially influences the solution 
of inflow region. 

One might think that the similarity of the 
boundary condition at ?/ = is not trivial, and 
the above boundary condition is imposed ad hoc. 
However, we must note that everything evolves 
self-similarly in this situation. In the self-similar 
stage VAot ^ D, the entire reconnection system 
including reconnection jet has no proper length 
other than the scale of FRWF Vao^. Thus, it 
naturally leads that the reconnection jet itself 
will grow self-similarly, and hence, the similarity 
growth of reconnection jet is realized. 

5.6. Relation to diffusion region 

The self-similar stage realizes after the system 
size VAot sufficiently exceeds any finite fixed spa- 
tial scale, i.e., the initial thickness of the current 
sheet and size of the diffusion region. Since both 
of these scales are estimated as D, the self-similar 
stage can be realized for Vao^ ^ -D as noted in sec- 
tion H There is no need to go into details about 
the structures of them because any fixed scale is 
negligible in the self-similar stage. In other words. 
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this shows the universahty of the self-similar evo- 
lution. The self-similar evolution does not de- 
pend on detailed structure of the current sheet. 
We can adopt any current sheet model if it satis- 
fies the initial MHD equilibrium. The self-similar 
evolution also does not depend on detailed resis- 
tivity model and structure of the diffusion region. 
The only function of the diffusion region which 
we need is to reconnect the field lines in a de- 
sired speed. We may expect that this is possible 
by dynamical change of the thickness of diffusion 
region. In fact, our numerical simulation showed 
that our new model is insensitive to the resistivity 
model (this is an important property of the fast 
reconnection). 

5.7. Comparison with previous steady and 
unsteady models 

There are several theoretical models for steady 
state and time-depending magnetic reconnection. 
We compare our self-similar evolution model with 
these previous models. Our discussion is focused 
only on fast reconnection, because very quick en- 
ergy conversion frequently observed in astrophys- 
ical phenomena suggests that fast reconnection 
should be considered as the responsible mecha- 
nism. 

The Petschek model (Petschek 1964) is charac- 
terized by a pair of slow shocks and a fast-mode 
rarefaction wave in the inflow region. This fast- 
mode rarefaction wave propagates outward from 
the reconnection point. As a result of the fast- 
mode rarefaction, the gradient of the magnetic 
field strength near the neutral point decreases due 
to the bending of magnetic field lines (Vasyliunas 
1975). This process suppresses and hmits the dif- 
fusion speed, and hence the reconnection rate. 

The Sonnerup model (Sonnerup 1970) is devel- 
oped from the Petschek model. This model is 
characterized by a pair of slow shocks and the hy- 
brid of fast-mode and slow-mode rarefaction. The 
fast-mode rarefaction wave is produced at the cen- 
tral region as in the Petschek model, but the slow- 
mode rarefaction wave is injected from the bound- 
ary. Because of the hybrid nature of rarefaction 
modes, the gradient of the magnetic field strength 
near the neutral point does not decrease as in the 
Petschek model. Therefore the Sonnerup model 
attains the maximum reconnection rate possible 
for magnetic energy converters (Priest & Forbes 
1986). A modified model of the Sonnerup solution 
was presented by Priest & Forbes (1986). They 



found a solution with a diffuse slow-mode rarefac- 
tion waves spread throughout the inflow region, 
while in the original Sonnerup model the rarefac- 
tion waves were treated as discontinuities. 

The property of the central region of our model 
can be categorized into the "almost-uniform re- 
connection" which means that magnetic field lines 
in the inflow region are almost straight and uni- 
form (see Priest & Forbes 2000). Priest & 
Forbes (1986) presented a unification scheme of 2- 
dimensional steady almost-uniform reconnection 
in a finite space filled with incompressible plasma. 
They found a family of continuous solutions char- 
acterized by a non-dimensional parameter "6o"- 
The Petschek {bo = 0) and Sonnerup (bo = 1) 
types are particular cases. Their result includes 
other types, e.g., Sweet-Parker type, fiux-pile-up 
type and stagnation-flow type according to the 
value of bo- These solutions {bo > 0) other than 
the Petschek type are characterized by the hybrid 
of slow-mode and fast-mode rarefaction. Our self- 
similar evolution is never influenced by the outer 
conditions, so the hybrid rarefaction never takes 
place in our case because there are no boundary- 
imposed waves in our case. We can conclude 
that the unique candidate which is worth com- 
paring with the central region of our model is the 
Petschek model. 

Figures |^ and ^ clearly show decreases of the 
magnetic field and gas pressure, thus the fast- 
mode rarefaction dominates in the vicinity of the 
diffusion region. Obviously the central region 
of this self-similarly evolving system is of the 
Petschek-type. Therefore, the reconnection rate 
might be limited in a way similar to the original 
Petschek model. This point will be clarified in our 
forthcoming paper. 

Far from the central region, there is a region in 
which the fast-mode compression takes place by 
the "piston effect" of the reconnection jet. This 
fast-mode compression causes the vortex-like re- 
turn fiow (see figure 13 of Nitta et al. 2001). The 
combination of the inner structure which resem- 
bles the original Petschek model and the outer 
structure involving the fiow vortices characterizes 
the evolutionary process, and the entire system 
unlimitedly expands self-similarly. 

Biernat, Heyn and Semenov (1987) and Se- 
menov et al. (1992) analytically studied the evo- 
lutionary process of fast magnetic reconnection 
in a similar situation to our case. Their analysis 
gives a general formalism of spontaneous time- 
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varying reconnection. They succeeded to ob- 
tain the solution for Petschek-type reconnection. 
However, in their works, plasma in the inflow re- 
gion is assumed to be incompressible for analyti- 
cal convenience. This means that the sound speed 
is infinitely large even when the Alfven speed is 
finite. This is equivalent to assuming that the in- 
flow region is filled with extremely high j3 plasmas 
(note that P ~ (sound speed) ^/ (Alfven speed) ■^). 
Needless to say, this assumption is unsuitable for 
most cases of astrophysical problems. Contrary to 
our case, fast-mode waves emitted from the cen- 
tral region instantly propagates to infinity. Al- 
though the proper spatial scale determined by 
the fast-mode wave propagation does not exist in 
their case, there is a possibility that the structure 
formed by slow-mode wave can be self similar, 
because the propagation speed of the slow-mode 
wave is finite in their case. Our work is under- 
stood as an extension of their works to realistic 
astrophysical systems filled with low /? compress- 
ible plasma. 

5.8. Astrophysical applications 

In the present level of our observational tech- 
nology, only one object in which we can obtain 
enough spatial and time resolution of the distri- 
butions of plasma parameters is solar fiares. We 
expect the application of our self-similar recon- 
nection model to the solar fiares. 

5.8.1. Dimming and inflow structure 

The propagation speed of FRW is estimated to 
be ~ 10^ [m s~^] for solar corona. Hence, the 
duration of the self-similar evolution is typically 
10^ — 10^ [s] for the fiux tubes having the spatial 
dimension ~ 10^ — 10^ [m]. The evolution hav- 
ing such a time scale is able to be resolved by the 
Solar-B project (required cadence of Solar-B X- 
ray telescope [XRT] is 2 sec, which is sufficient 
to resolve the expected time evolution of the self- 
similar reconnection, see Golub 2000.). Especially 
we expect that the "dimming" will be detected by 
Solar-B. Expected X-ray image of the dimming 
is argued in Nitta et al. 2001 (see section 4.5). 
We consider that the dimming is a naturally ex- 
pected phenomenon at the infiow region of our 
self-similar reconnection model. If we obtain the 
information of velocity field and other plasma pa- 
rameters at the reconnection site by observation, 
we can inspect our model more directly. 



The Solar-B satellite will clarify also the veloc- 
ity field, plasma density and temperature near 
the reconnection point by the extreme ultravio- 
let imaging spectrometer (EIS). The combination 
of XRT and EIS will reveal a detailed structure of 
the reconnection system in solar fiares. Yokoyama 
et al. (2001) found the evidence for reconnection 
inflow in a flare on 1999 March 18. The estimated 
inflow speed is 5 km s^^ which corresponds to the 
reconnection rate 0.001-0.03 (this result is not far 
from the expected value 0.05 from our model). 
Their result encourages our theoretical effort of 
applying our self-similar model to the infiow re- 
gion of real solar fiares. We hope that the ex- 
pected feature of the infiow region by virtue of our 
self-similar reconnection model will be inspected 
by the Solar-B project. 

5.8.2. Application to realistic cases 

Roussev et al. (2001 a, b, c) predicted emis- 
sion from explosive events in the solar transition 
region induced by magnetic reconnection. Tlicy 
performed MHD numerical simulations including 
thermal conduction, radiative losses and volumet- 
ric heating. They treated spontaneous evolution 
of magnetic reconnection with the initial condi- 
tions similar to our model (anti-parallel magnetic 
field and asymptotically uniform plasma distribu- 
tion) and with the resistivity (artificially localized 
and kept as a constant) that is also very similar 
to ours. The resultant dynamical properties are 
similar to our study, indicating that the essence of 
our model is applicable to more complicated sit- 
uations with more realistic thermal and radiative 
properties. 

5.9. Expectation of 3D self-similar reconnection 

Our present discussion is restricted to 2D re- 
connection in this as well as in the previous paper 
(Nitta et al. 2001) for simplicity. However, actual 
current sheet systems have finite depth in the di- 
rection perpendicular to the figures of this paper 
(for example, radius of the fiux tube or length of 
the arcade structure of a bunch of fiux tubes). 
We believe the essence of the self-similar growth 
has already been understood by our 2D approx- 
imation. When we discuss the phenomenology 
of the evolutionary reconnection by virtue of the 
self-similar reconnection model in near future, we 
need a more precise 3D model. 

We can expect that the self-similar growth also 
takes place in a 3D system. As discussed in sec- 
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tion |5]^, a sufficiently evolved system has only 
one proper spatial scale, i.e., the scale of FRWF 
VAot which is increasing in proportion to time 
(note that the propagation speed of the fast-mode 
is almost isotropic in a 3D current sheet system 
filled with a low (3 plasma, hence the shape of 
FRWF will be spherical). This is just the situ- 
ation that leads to realize the self-similar expan- 



sion. Of course, detailed structure of the 3D re- 
connection system will be different from our 2D 
model, but essential properties will be common 
between them. 

One of the author SN is grateful to T. Kudoh at 
National Astronomical Observatory of Japan for 
useful comments and fruitful discussion. 



APPENDIX 

APPENDIX A: DERIVATION OF MHD EQUATIONS IN THE ZOOM-OUT COORDINATE 

We show the derivation of equations briefly in this appendix. The zoom-out coordinate is 

deflned by equation ([1|). The velocity in this new coordinate has been given in equation (^. This 
equation is deformed to 

Vaov' = v- Vaot' . (Al) 

The second term of the right hand side is the apparent converging flow properly appearing in the 
zoom-out coordinate. Physical quantities are put into non-dimensional forms as in (D-(D. With these 
preparations, we can derive the modifled MHD equations in the zoom-out coordinate. 

Any scalar fleld, e.g., f{r,t) is transformed in the zoom-out coordinate as f{r'[r,t],t). Thus the time 
derivative of any scalar fleld can be transformed as 



dl 

dt 



df 




Or' 


dt 


r' 




df 






dt 


r' 


t 



V'f (A2) 
V7 (A3) 



where subscripts r and r' show the time derivative in the flxed coordinate and zoom-out coordinate, 
respectively, and V' denotes the spatial derivative with respect to r'. Here we have used (0) to derive 
the second equality. The transform of time derivative for any vector fleld is similar to this. 
The gradient of any scalar fleld is transformed as 

V/lr- = :^V7 . (A4) 

The divergence and rotation of any vector fleld, e.g., A are transformed as 

V-A^J-V-.A (A5) 

and 

V X A = X A , (A6) 

respectively. By using these transformed notations, we can rewrite the MHD equations as follows. 

dp 

^ + V ■ (pv) = (continuity eq.) (A7) 

+ V ■ (pvv - —Bb] + V (p + —\ = (momentun eq.) (A8) 
dt \ Ho I \ 2/io/ 
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dB 

— + V ■ {vB - Bv) = (induction eq.) (A9) 

Pp^'^ = Const, (polytrope relation) (AlO) 

As a representative case, we will demonstrate the derivation of the modified continuity equation. The 
derivations for other modified MHD equations are quite similar to this. 
The original continuity equation 

dp 
dt 

is transformed to 



+ V ■ (pv) = (All) 



r 



dp 
dt 



- y • V'p + T^V' ■ (pv) = . (A12) 



Here we suppose the stationariness in the zoom-out coordinate (self-similar expansion). It makes the 
first term vanish. By substituting the normalization relations (^-(^D into this transformed equation, 
we easily obtain the following equation, 

r' ■ Vy - V ■ [p'{v' + r')] = . (A13) 

This equation is expanded to the form 

- vy ■ v' - p'V' • v' - p'V ■ r' = . (A14) 

We combine the first two terms to obtain —V' ■ {p'v'). We should note that in the two-dimensional 
problem, V' ■ r' = 2. Thus we finally obtain the transformed continuity equation (^. 

Similarly we can easily obtain the other modified equations (0)-(0)- Each derivation is not difficult, 
but is very tedious, hence we omit detailed derivation here. 

APPENDIX B: PERTURB ATIVE RELATIONS 

By assuming that the deviation by reconnection from the initial equilibrium state is very small, we 
derive the linearized relations. We approximately expand each quantity as in equations ([I^)-(p!7|) with 
(p!0D-(p!3D, and substitute them into the modified MHD equations (|)-(|)- We can obtain the linearized 
MHD equations as is demonstrated below. 

For example, we will demonstrate the derivation of the linearized continuity equation (p^). The 
derivation of other linearized MHD equations is quite similar. We substitute the expansion forms (|l^ 
and (pTsD into the modified continuity equation (|). By expanding equation (|^) up to the first order of 
magnitude, we obtain two equations: the zeroth order equation 

V ■ (p>o') = -2p'o (Bl) 

and the first order equation 

V ■ {p',v^' + p[vo') = -2p[ . (B2) 
The left hand side of the zeroth order equation is reduced to 

V ■ (p>o') = V ■ i-r') (B3) 
= -2 . (B4) 

The right hand side of this equation is -2. Thus, we can see that the zeroth order equation is trivial and 
has no information. The left hand side of the first order equation is reduced to 

V ■ (p>i' + p>o') = V ■ V,' + V ■ {-p[r') (B5) 
= V ■ Vi' - r' ■ vy - p'lV • r' . (B6) 

By noting that the last term is — 2p'^ and together with the right hand side, we obtain the linearized 
continuity equation (|18|). By similar reduction, we obtain other linearized equations (p!9D-(pl|). 



Nitta et al. 



13 



APPENDIX C: DERIVATION OF THE G-S EQUATION AND OTHER RELATIONS 

We here show the reduction of (^31)- (^61). Let us start with the deformation of the hnearized momentum 
equation (P^. Each term of that equation is reduced as following: 



- i ■ V'B[ 



V'{i ■ B[) 



dxdy dx'^ 



r ■ V 



X- 



Ix 



dx 



dxdy 
dy 



dy^ 



y- 



X- 



dx 



+ y- 



dv 



dy 



Then we obtain the reduced form of each component of the equation (|T 

dv'. 



as 



X- 



dx 



+ y- 



dv'. 



X- 



1?/ 



dx 



+ y- 



dv\ 



dy 



+ 



dy 
( d^ 



0, 



+ 



d^ 



dx'^ dy 



A\ =0 



Similarly, from the linearized induction equation (|20|), we obtain the following two equations: 



d^A' d^A\ 



dv'^ 



x- 



dxdy 
d^A[ 



+ y 



dy"^ 
d^A'^ 



dy 



dv'. 



. 



(CI) 

(C2) 
(C3) 

(C4) 
(C5) 

(C6) 
(C7) 



dx"^ ' " dydx dx 

We should note that these two equations from the induction equation have the same information in 



common 



+ y^ - K - <y 



dx 



dy 



(Const, independent of x, y) 



(CJ 



The vector potential A'^ has a freedom of an indeterminate constant, then the constant on the right 
hand side can be absorbed into the potential, and we obtain equation (^If). 

We substitute the two components of the induction equation (|C6D and (|C7|) into (|C5|) , and obtain a 
single partial differential equation (^) for A\ after simple calculations. This is the G-S equation. 

Another component of the momentum equation (|C4|) is transformed to a simple form: 

r' ■ VS. = . (C9) 

This equation clearly shows that v'^^ depends only on the angle 6 from the x-axis. We should remember 
that v'l^ = on r' = 1 (FRWF). Thus this equation implies v'^^ = everywhere in our linearized analysis. 
This is the equation (^). 

The linearized continuity equation (|18D is deformed by using the above results. The first term is 
reduced by the help of ( |C6|) as 



dv[ 



x 



dy 
d^A\ 



+ y 



d^A'^ 



dxdy dy 



r'. V 



dy 



(CIO) 
(CU) 
(C12) 
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Together with the second term of (|T8D, we can see that the continuity equation implies 

^'■V'(^-p;)=0. (C13) 

Similarly to the discussion of the previous paragraph, we obtain the result ( pS]) because we have set the 
boundary condition dA[/dy = and p[ = on FRWF. The relation ( p6D is easily obtained from (|2T|). 

APPENDIX 
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Fig. 1. — Boundary condition of A[ along the slow shock {y = 0). The dotted line shows the result of our simulation 
in Nitta et al. 2001 (Case A). The solid line is the approximated boundary condition adopted in this work. We can see 
the approximated model is very simple and gives good agreement with the simulation result. 

Fig. 2. — Contours of the perturbed flux function (z-component of magnetic vector potential) A[. The solid (dashed) 
contours correspond to positive (negative) value of A[ (The same applies to figures ^, ||, ||, || and |^). The contact 
discontinuity is located at a; = Xc ~ 0.64, y = 0. The location of the maximum value is at the reconnection point 
(x = y = 0), and that of minimum value is near the spearhead of reconnection jet {x = Xm ^ 0.84, y = 0). Hence the 
first order magnetic field is turned counter-clockwise along the left-half arcs, and clockwise in the right-half arcs. The 
zeroth order magnetic field is directed to the right in the upper-half plane. Thus the strength of the magnetic field is 
reduced by the first order variation in the vicinity of the reconnection point. This region is filled with the fast- mode 
rarefaction dominated inflow. This property is common to the original Petschek model. On the contrary, the strength 
of the magnetic field is enhanced by the first order variation near the spearhead of reconnection jet. This is caused by 
the fast-mode compression (piston effect) produced by the reconnection jet. This compressed region is characteristic to 
our dynamically evolving system, and quite different from previous stationary models. 

Fig. 3. — The variation of flux function from the initial state (result of our simulation in Nitta et al. [2001]). We 
can see essentially the same feature in figures ^ and |^ except that the configuration near the cc-axis (y < 0.05) is rather 
different. This is due to the existence of reconnection jets which are neglected in our present analytical study. 

Fig. 4. — Contours of the y-component v'ly of the velocity perturbation. The a;-component of the velocity perturbation 
is zero in the linear analysis. We should note that the value of v[y is equivalent to the strength of electric field in the z 
direction (so-called reconnection electric field) in the normalized unit. The value near the reconnection point is negative. 
This means that the inflow is sucked into the reconnection jet. Positive values exist near the right-half outer boundary. 
This is the outflow gushing out from the region near the spearhead of the reconnection jet. Hence, this outflow is a 
secondary flow properly characteristic to the evolutionary process. This outflow enhances the vortex-like flow. The value 
at the reconnection point {x = y = 0) represents the reconnection rate of the system, and the value is 0.055 from our 
simulation (Nitta et al. 2001). 

Fig. 5. — The inflow velocity (y-component) (result of our simulation if Nitta et al. [2001]). 
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Fig. 6. — Contours of the perturbed density p\. If we divide this value by adiabatic constant 7, these contours 
show the perturbation in gas pressure P'^ in the normahzed unit. The density (hence pressure) is reduced near the 
reconnection point, and enhanced near the spearhead of the reconnection jet. These features are also caused by the 
fast-mode rarefaction and compression as discussed in the caption for figure ^. 

Fig. 7. — The variation of density distribution from the initial state (result of our simulation in Nitta et al. [2001]). 
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